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SUMMARY 


This  paper  treats  the  effect  of  wind-tunnel  walls  on  the  oscillating 
two-dimensional  air  forces  in  a  compressible  medium.  The  walls  are  simu¬ 
lated  by  the  usual  method  of  placing  images  at  appropriate  distances 
above  and  below  the  wing.  An  important  result  shown  is  that,  for  certain 
conditions  of  wing  frequency,  tunnel  height,  and  Mach  number,  the  tunnel 
and  wing  may  form  a  resonant  system  so  that  the  forces  on  the  wing  are 
greatly  changed  from  the  condition  of  no  tunnel  walls.  It  is  pointed 
out  that  similar  conditions  exist  for  three-dimensional  flow  in  circular 
and  rectangular  tunnels  and  apparently,  within  certain  Mach  number  ranges, 
in  tunnels  of  nonuniform  cross  section  or  even  in  open  tunnels  or  jets. 


INTRODUCTION 


The  understanding  of  flutter  and  other  nonsteady  phenomena  requires 
a  knowledge  of  the  associated  unsteady  flow.  In  the  underlying  theories 
of  unsteady  flow,  such  assumptions  as  small  displacements,  linearizations, 
and  an  inviscid  fluid  are  made  in  order  to  obtain  workable  and  usable 
results.  When  it  is  necessary  to  investigate  the  effect  of  these  assump¬ 
tions  on  analytical  results  by  measurements  of  the  forces  and  moments  on 
an  oscillating  wing  in  a  wind  tunnel  or  to  treat  cases  that  do  not  con¬ 
form  to  theory,  the  question  of  the  effect  of  the  tunnel  walls  naturally 
arises .  In  the  case  of  steady  flow  the  problem  of  the  effect  of  tunnel 
walls  is  more  or  less  classic  and  has  been  treated  by  many  investigators. 
In  general,  these  investigators  have  been  able  to  obtain  relatively  simple 
factors  which  can  be  used  to  modify  measurements  of  the  air  forces  on  a 
wing  in  a  tunnel  to  correspond  to  free-air  conditions.  The  extension 
of  the  results  to  compressible  flow  presents  no  difficulties  since  the 
results  for  incompressible  flow  can  be  corrected  according  to  Prandtl- 
Glauert  correction  factors. 
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In  the  case  of  unsteady  flow,  Reissner,  reference  1,  and  W.  P.  Jones, 
reference  2,  have  published  papers  showing  the  effect  of  wind-tunnel  walls 
for  the  incompressible  case.  In  both  papers,  the  influence  of  the  tunnel 
walls  is  found  to  be  comparatively  small  for  most  cases,  although  indica¬ 
tions  are  given  that,  for  some  ranges  of  a  reduced-frequency  parameter, 
the  effect  may  be  quite  large.  In  the  unsteady  case,  unlike  the  steady 
case,  the  transition  from  results  for  incompressible  flow  to  those  for 
compressible  flow  cannot  be  accomplished  by  simple  transformations.  This 
difficulty  is  a  result  of  the  fact  that,  in  an  incompressible  fluid,  the 
velocity  of  propagation  of  a  disturbance  is  infinite  and  no  time  lag 
occurs  between  the  initiation  of  a  disturbance  and  its  effect  at  another 
position  in  the  field,  but,  in  a  compressible  fluid,  a  definite  time  is 
required  for  a  signal  to  reach  a  distant  field  point  so  that  both  a  phase 
lag  and  a  change  in  magnitude  result.  Under  certain  conditions  this 
phase  lag  can  result  in  a  resonant  condition  which  would  involve  large 
corrections . 

The  purpose  of  this  paper  is  to  consider  the  effect  of  wind-tunnel 
walls  on  the  forces  on  an  oscillating  airfoil  of  infinite  span  with  con¬ 
siderations  of  the  compressibility  of  the  fluid.  The  usual  method  of 
images  is  employed  in  order  to  satisfy  the  condition  of  no  normal  velocity 
at  the  tunnel  walls.  First,  the  effect  of  tunnel  walls  on  the  induced 
vertical  velocity,  hereinafter  referred  to  as  downwash,  of  an  oscillating 
doublet  is  determined  and  this  result  is  used  to  formulate  the  integral 
equation  for  the  downwash  of  an  oscillating  airfoil  in  a  tunnel.  This 
paper  is  not  intended  to  give  numerical  values  or  any  detailed  calcula¬ 
tions  of  final  tunnel-wall  correction  factors  but  mainly  to  show  the 
existing  need  for  such  calculations  and  to  present  equations  for  calcu¬ 
lating  corrections  for  the  two-dimensional  case. 


SYMBOLS 


A 

b 

c 

H 

Hq(2),  %(*> 

M 


constant 
semichord 
velocity  of  sound 
tunnel  height 
Hankel  functions 
Mach  number 


Ap 


local  pressure  difference 
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tijne 


V 

w,  w0,  w-l 
t,  x,  y,  z, 
P  =  /l  -  M2 

y 

CO 

X 

t 

9  . 

P 


velocity 

downwash  or  vertical  induced  velocity 
Cartesian  coordinates 

Euler*  s  constant 
angular  frequency 
wave  length 
acceleration  potential 
velocity  potential 
fluid  density 


ANALYSIS 

Effect  of  Tunnel  Walls  on  the  Downwash  of  a  Single  Doublet 


The  differential  equation  that  governs  flow  due  to  small  nonsteady 
perturbations  imposed  on  a  steady,  uniform  flow  field  is  the  wave  equa¬ 
tion.  Referred  to  rectangular  coordinates,  fixed  relative  to  the  undis¬ 
turbed  stream  at  infinity,  this  equation  is 


(1  -  M2)  + 

dx2 


2M  a2^ 

C  dx  dt 


1  ii.o 

c2  St2 


(1) 


In  this  equation  the  independent  variable  \|r  may  be  regarded  as  either 
a  perturbation  velocity  potential  or  as  an  acceleration  potential.  In 
treating  the  boundary  conditions  of  the  second  section  of  this  analysis 
it  is  convenient  to  regard  y  as  an  acceleration  potential.  Thus,  in 
order  to  be  consistent,  t  is  hereinafter  regarded  as  an  acceleration 
potential.  Accordingly  i|r  is  directly  proportional  to  a  perturbation 
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pressure  field  and  is  therefore  related  to  a  perturbation  velocity 
potential  cp  as  follows: 


f  = 


3®  3cp 

— ■  +  V  — 

3t  3x 


(2) 


In  order  to  calculate  the  downwash  w  =  —  associated  with  it  is 

3y 

necessary  to  solve  equation  (2)  for  cp  in  terms  of 

When  f  and  cp  are  sinusoidal  functions  of  time,  such  that 


^x,y,t)  =  e^^Xjy) 
<P(x,y,t)  =  e1“t9(xJy)J>" 


(3) 


equation  (2)  becomes  independent  of  time  and  thus  reduces  to  an  equation 
with  one  dependent  variable *  namely 


=  iaxp  +  V 


8cp 

dx 


(h) 


This  equation  can  be  integrated  with  respect  to  x  to  give 


(5) 


where  the  lower  limit  of  integration  is  chosen  for  later  convenience  so 
that  cp  vanishes  far  ahead  of  the  point  of  disturbance.  The  downwash 
may  be  readily  calculated  with  the  use  of  this  equation.  In  the  absence 
of  tunnel  walls  the  retarded  potential  \|^q  (that  is,  the  potential  corre¬ 
sponding  to  outgoing  waves)  of  a  harmonically  pulsating  pressure  doublet 
located,  for  simplicity,  at  (0,0)  that  satisfies  equation  (l)  is 
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where  Hg^  and  are  Hankel  functions  as  defined  in  reference  3, 

A  is  an  arbitrary  constant  denoting  doublet  strength }  oo  is  circular 

frequency,  and  p  =  |/l  -  M2.  The  Hankel  function  %(?)  in  equation  (6) 
becomes  infinite  /as  — ^  a>s  its  argument  approaches  zero.  Other- 


w^-se  ^1  ^-S  continuous  and  approaches  zero  as  its  argument  approaches 

infinity.  Thus  the  only  discontinuity  in  is  at  the  location  of  the 

doublet,  that  is,  at  (x=0,y=0). 


In  the  presence  of  plane  tunnel  walls  located  parallel  to  the  x-axis 
at  H/2  units  above  and  H/2  units  below  the  doublet  position,  the 
potential .  f  of  a  pressure  doublet  may  be  represented  by  the  potential 
of  an  infinite  system  of  appropriately  chosen  reflecting  doublets,  namely 
(see  fig.  1)  3 


In  this  equation  the  term  corresponding  to  n  =  0  is  the  potential 

equation  (6),  discussed  in  the  preceding  paragraph.  It  may  be  noted  that 
only  this  term  of  the  infinite  summation  in  equation  (7)  gives  rise  to  a 


discontinuity  in  f  at  any  point  within  the  tunnel  |^y^| <x<o^. 
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The  infinity  of  terms  corresponding  to  n  f  0  are  necessary  to  cause 
the  downwash  w  to  vanish  at  all  points  of  the  tunnel  walls,  y  =  ±|. 

The  downwash  along  the  midsection  of  the  tunnel  y  =  0  is  given  by 


ico£ 

(-i)V?  4  “0(2) 

n=-oo  8y 


6) 


cp: 


]p  +  P2(y  -  nH)2 


d| 


=  w0  +  wx 


(8) 


where 


represents  the  downwash  associated  with  the  pressure  doublet  in  the 
absence  of  tunnel  walls  and 


ico| 

C-1)^  4  -0<2> 

sy2 


CO 


cp 


/l2  +  P2(y-nH): 


d| 


(10) 

represents  the  additional  downwash  due  to  the  presence  of  tunnel  walls. 
Thus  the  relative  value  of  Wq  as  compared  with  wq  +  wq  is  the  main 
item  of  interest  here. 

The  integrals  appearing  in  equations  (9)  and  (10)  can  be  reduced  to 
simpler  form  for  evaluation  but  since  the  steps  required  to  reduce  one 
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of  the  integrals  are  the  same  as  required  to  reduce  the  other,  only 
the  integral  appearing  in  equation  (9)  will  be  treated  in  detail.  The 
reduced  form  of  the  other  integral  can  then  be  obtained  by  simple  com¬ 
parison.  The  Hankel  function  in  equation  (9)  satisfies  the  following 
identity: 


Substituting  this  relation  into  equation  (9)  gives 


In  equation  (12)  the  first  integral  can  be  integrated  twice  by  parts  to 
give  for  wq 


(13) 
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By  writing  the  integral  in  equation  (13)  as  the  sum  of  two  integrals, 
namely 


+ 


and  making  a  change  of  variable 


the  expression  for  Wq  may  be  further  reduced  to 


(111) 


(15) 


In  the  limit  y  =  0  the  expression  in  braces  in  equation  (16)  reduces 
to  the  kernel  of  Possio's  integral  equation  relating  pressure  and  down- 
wash  for  the  oscillating  airfoil  in  compressible  flow.  (This  result 
checks  the  results  for  this  expression  given,  for  example,  in  reference  H.) 
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In  general^  this  infinite-series  representation  of  w-j_,  equation  (17), 
converges  to  a  finite  value.  However;  for  certain  critical  values  of  the 
frequency  parameter  oaH/V,  it  is  found  that  the  value  of  w^  becomes 
infinite.  This  fact  can  be  readily  made  evident  by  use  of  relations 
given  in  reference  5  where  it  is  shown  that  an  infinite  series  of  Hankel 
functions  of  the  type  appearing  in  equation  (17)  can  be  replaced  by  an 
equivalent  series  of  exponential  functions  as  follows: 


00 

-  —I 

00 

/  ✓  .2 

JL  (-dV2> 

-%/x2  *  (fnH)2 

=  Y_  (-l)nH0(2) 

*  (i) 

n=l 

Lsr 

n=l 

Vp  f  \pH/ 

It  may  be  seen  that,  if  this  relation  is  substituted  into  equation  (17), 
the  value  of  w^  becomes  infinite  for  all  values  of  x  when  the  fre¬ 
quency  parameter  coH/V  has  any  of  the  values  given  by 


wH 

V 


-  (2m  -  1)  ^ 


(m  -  1,  2,  3,  . 


(19) 
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These  critical  values  of  the  frequency  parameter  correspond  to  a  condi¬ 
tion  of  pure  resonance  in  the  tunnel  which  in  the  present  case  implies 
that  a  harmonic  disturbance  of  any  finite  amplitude  may  lead  to  a  down- 
wash  of  infinite  amplitude. 

Of  course  these  infinite  values  of  w-^  would  never  be  realized 
under  practicable  conditions  because  factors  such  as  finite  tunnel  length, 
absorption  through  walls,  fluid  viscosity,  and  so  forth  that  would  give 
rise  to  damping  would  make  pure  resonance  unobtainable;  however,  with 
damping  present,  resonant  frequencies  yielding  values  of  ccH/V  in  the 
neighborhood  of  those  given  in  equation  (19 )  would  exist  and  it  is  not 
likely  that  quantitative  agreement  or  even  possibly  qualitative  agreement 
between  calculated  and  measured  downwash  (or  forces)  can  be  realized  when 
the  value  of  coH/V  is  in  the  neighborhood  of  these  critical  values. 

It  is  interesting  to  note  that  the  effect  of  boundary  conditions 
such  as  section  geometry,  tunnel-wall  flexibility,  and  so  forth  is  to 
change  the  value  of  the  critical  frequency  but  not  to  do  away  with  the 
possibility  of  resonance.  Also,  by  treatments  similar  to  those  employed 
herein,  it  can  be  shown  that  under  idealized  conditions  resonance  can 
occur  in  three-dimensional  flow  in  both  round  and  rectangular  tunnels  or 
apparently,  within  certain  Mach  number  ranges,  in  tunnels  of  nonuniform 
cross  section  (expanding  or  contracting  section)  or  even  in  open  tunnels 
or  jets. 

The  fundamental  or  smallest  critical  values  of  ccH/V,  corresponding 
to  m  =  1  in  equation  (19),  are  shown  plotted  as  functions  of  Mach  num¬ 
ber  M  in  figure  2.  This  figure  indicates  that  there  is  no  finite 
critical  value  of  coH/V  for  the  conditions  M  =  0,  V  ^  0,  and  c  =  00, 
which  correspond  to  a  flow  of  incompressible  fluid  in  the  tunnel.  This 
result  agrees  with  those  found  in  references  1  and  2. 

The  frequency  parameter 


^  =  (2m  -  l)np  (m  =  1,  2,  3,  .  .  .)  (20) 

c 


which  may  be  derived  from  equation  (19)  is  shown  plotted,  for  m  =  1, 
as  a  function  of  Mach  number  in  figure  3-  Equations  (20)  and  figure  3 
show  that  finite  values  of  the  critical  frequency  exist  for  the  condi¬ 
tions  M  =  0,  V  =  0,  and  c  £  00 .  These  conditions  correspond  to  a  com¬ 
pressible  fluid  at  zero  velocity  in  the  tunnel.  For  these  conditions 
equations  (20)  and  the  corresponding  wave  lengths 


2H 


.)  (21) 


2m  -  1 


(m  1,  2,  3j  •  • 
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agree,  respectively,  with  results  found  in  the  literature  for  the  char¬ 
acteristic  frequencies  and  wave  lengths  associated  with  transverse 
acoustic  vibrations  in  rectangular  chambers  when  the  location  of  the 
source  of  disturbance  is  excluded  as  a  nodal  point.  See,  for  example, 
reference  6. 

It  may  be  of  interest  to  note  that  equation  (20)  can  be  derived 
from  the  principle  of  standing  waves  as  follows:  The  condition  for 
resonance  for  the  type  of  disturbance  considered  implies  that  the 
standing  transverse  waves  have  a  maximum  velocity  at  the  midsection  of 
the  tunnel  and  zero  velocity  at  the  boundaries.  A  half-sine  wave  of 
wave  length  X  =  2H  or  any  odd  divisor  of  this  length,  namely, 

2H 

X  =  - - -  satisfies  this  condition.  If  c  is  the  velocity  of  sound 

2m  -  1 

in  the  medium  and  V  the  velocity  of  the  medium,  the  velocity  of  pro- 
pagation  of  a  disturbance  in  a  fixed  plane  perpendicular  to  the  air  flow 

is  .  Since  the  frequency  is  given  by  the  speed  of  propagation 

divided  by  the  wave  length  there  is  obtained 


f  = 


(2m  -  l)/c2  -  V2 
2H 


CO 

2tt 


or 


f  =  up (2m  -  1) 


Integral  Equation  for  an  Airfoil  of  Infinite  Aspect 
Ratio  Oscillating  in  a  Wind  Tunnel 

In  order  to  present  equations  from  which  tunnel-wall  corrections 
for  two-dimensional  flow  can  be  calculated,  use  is  made  of  the  foregoing 
analysis  to  derive  the  integral  equation,  relating  downwash  distributions 
and  lift  distributions,  for  the  effect  of  tunnel  walls  on  the  lift  dis¬ 
tribution  associated  with  a  given  downwash  distribution. 

The  resultant  pressure  or  local  lift  Ap  associated  with  the  accele¬ 
ration  potential  of  a  single  doublet  located  at  (xq,0)  with  strength 

depending  on  streamwise  position  Xq  may  be  expressed  simply  as  (compare 
with  equation  (6)): 
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For  an  airfoil  inside  a  two-dimensional  tunnel  the  relation  between  downwash  and  local  lift 
becomes  (compare  with  equation  (8)) 
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results  can  be  obtained  for  these  cases 
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CONCLUDING  REMARKS 


The  important  result  shown  is  that,  in  a  tunnel  of  infinite  length 
containing  a  flowing  fluid,  a  resonant  condition  involving  a  transverse 
oscillation  of  the  fluid  across  the  tunnel  is  possible  and  measured  air 
forces  at  or  near  this  condition  of  resonance  might  be  greatly  modified 
from  those  measured  in  free  air.  This  resonant  condition  is  a  (simple) 
function  of  Mach  number,  tunnel  height,  and  wing  frequency  and  brings  to 
attention  a  new  type  of  tunnel-wall  interference. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  September  2k,  1951 


REFERENCES 


1.  Reissner,  Erics  Wind  Tunnel  Corrections  for  the  Two-Dimensional 

Theory  of  Oscillating  Airfoils.  Rep.  No.  SB-3I8-S-3,  Cornell  Aero. 
Lab.,  Inc.,  April  22,  I9I4.7. 

2.  Jones,  W.  Prichard:  Wind  Tunnel  Interference  Effect  on  the  Values  of 

Experimentally  Determined  Derivative  Coefficients  for  Oscillating 
Airfoils.  R.  &  M.  No.  1912,  British  A. R.C. ,  19^3. 

3.  Watson,  G.  N.s  A  Treatise  on  the  Theory  of  Bessel  Functions.  Second 

ed..  The  Macmillan  Co.,  I9I1U. 

U.  Karp,  S.  N.,  Shu,  S.  S.,  and  Weil,  H.s  Aerodynamics  of  the  Oscillating 
Airfoil  in  Compressible  Flow.  Tech.  Rep.  No.  F-TR-II67-ND,  Air 
Materiel  Command,  U.  S.  Air  Force,  Oct.  19k7 . 

5.  Inf  eld,  L.,  Smith,  V.  G.,  and  Chien,  W.  Z.s  On  Some  Series  of  ne -.sel 

Functions.  Jour.  Math,  and  Phys.,  vol.  XXVI,  no.  1,  April  19kl . 

6.  Morse,  Philip  M.s  Vibration  and  Sound.  Second  ed.,  McGraw-Hill  Book 

Co.,  Inc.,  I9U8. 


1 6 


MCA  TN  2552 


Figure  1.-  Sketch  showing  reflecting  system  of  doublets  simulating 

two-dimensional  tunnel  walls. 
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Figure  3.-  Fundamental  critical  values  of  frequency  parameter  aE/c 
plotted  as  a  function  of  Mach  number  M. 
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